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The Kwiecihski evolution of unintegrated parton distributions (UPDs) in the 
transverse-coordinate space is analyzed with the help of the Mellin transform. Nu- 
merical results are presented for the unintegrated pion distributions with a simple 
valence-like initial condition at the low scale, which follows from chiral large- N c 
quark models. The effect of spreading of UPDs in the transverse momentum with 
the increasing scale is confirmed, with (kj_) growing asymptotically as Q 2 a(Q 2 ). 
Formal aspects of the equations, such as the limits of UPDs at x — ► 0, x — ► 1, and 
at low and large 6, are straightforward to obtain with our method. 

This talk is based on Refs. [1,2] and focuses on numerical as weel as formal 
aspects of the QCD evolution of UPDs proposed by Kwiecihski [3]. Practical ap- 
plications are presented by Szczurek [4], while similar approaches are discussed by 
Jung [5] and Lonnblad [6] in these proceedings. 

The definition of the leading- twist UPDs [7] may be read-off from Fig. 1(a). 
The transverse momentum in the partonic loop, k±, is left unintegrated and is not 
constrained. Kwieciriski's approach is based on the CCFM formalism with the 
following modifications: 

• One-loop CCFM is used, with the angular ordering replaced with a stronger 
condition on the scaled gluon momenta of the form q' ±i > ^_ l5 where 
q[ = qij(\ — Zi), Zi = Xi/xi-i, and Xi denotes the fraction of the longitu- 
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Figure 1. (a) Definition of the leading- twist UPD. (b) Kinematics inside the partonic cascade. 
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dinal momentum of the hadron carried by the parton. 

• Both singlet and non-singlet quarks are included. 

• Non-Sudakov form factor is set to unity. 

• The evolution is written in the transverse-coordinate space, 6, Fourier- 
conjugated to the transverse momentum. 

The Kwiecihski equations are diagonal in the 6-space and read [3]: 

W dQ 2 " 1 } 



g s (Q 2 ) 

2tt 



2\ /•! 



J dzP qq (z) e(z-x)J ((l-z)Qb)f NS (^b,Q^-hs(x,b,Q) 



with similar coupled equations for the sea quarks and gluons. Note that for the 
integrated PDs (case of b = 0) the equations assume the familiar DGLAP form. 

In this talk we consider the uPDF of the pion. The initial conditions for the 
evolution are taken from the chiral quark models [2,8,9], which in the chiral limit 
provide a particularly simple, factorized form. Moreover, these conditions are con- 
sistent with the formal requirements, such as proper normalization, correct support, 
crossing symmetry, and the momentum sum rule. The working scale of the model, 
Qo, is estimated using the fact that at Qo all the momementum is carried by the 
quarks, which are the only degrees of freedom at that scale (see [2] for details). The 
Spectral Quark Model [9] provides a factorized initial condition of the form 

q(x,b,Q o )=F(b)0(x)0(l-x), (2) 

M v b \ 
2 J 



F(6)=(l + ^)ex P 



for the valence quarks, with no gluons nor sea quarks at the scale Qo- Note that the 
functional dependence in b in Eq. J2J) factorizes in Eq. QJ. Thus, it is convenient to 
denote f(x, 6, Q) = F(b)f evol (x, 6, Q) and present the results for the the evolution- 
generated part / evo1 . 

Equations can be solved very efficiently using the Mellin transform [2]. In 
the moment space we find 

, 2 #N s(n,fr,Q) = a(Q 2 ) 
dQ 2 8tt 

where f^s(n,b,Q) denote the nth moment, and 7 n ,A/ r s'(fr) are the 6-dependent 
anomalous dimensions [2] . The formal solution of Eq. (|3J) is 

f Q2 dQ' 2 a(Q' 2 ) ~ 



Q ^ = ^^In^siQbjfNsin, b, Q), (3) 



/ NS (n,6,Q) 

— exp 



fas(n,b,Q ) 



(4) 



and similarly in the singlet channel. In ordert to come back to the x-space, the 
inverse Mellin transform is carried out numerically. The numerical results are shown 
in Figs.OandEl In Fig.^we plot the mean squared transverse momentum generated 
by the evolution. As Q is increased, (k\) grows as ^(Q 2 )^ 2 . 
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Figure 2. Evolution-generated UPDs' for the°pion for various values &f trie transverse coordinate 
(from top to bottom: 6 = 0, 1, 2, 3, 4, 5 and 10 GeV -1 ), plotted as functions of the Bjorken x. 
The evolution is made with the initial condition J2J at Qo = 313 MeV up to Q = 2 GeV. Solid 
lines - gluons, dashed lines - valence quarks, dotted lines - sea quarks. 
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Figure 3. Evolution-generated UPDs for the pion for Q = 2 Gety[Qf$~x]= 0.1, plotted as functions 
of 6. The evolution is made with the initial condition J2j at Qo = 313 MeV. The numerical results 
are represented by squares for the non-singlet quarks, diamonds for the singlet quarks, and stars for 
the gluons, while the solid line shows the asymptotic power-law formula for the case of non-singlet 
quarks. We note the much faster fall-of for the gluons than for the quarks. As Q is increased or 
x decreased, the distributions in b become narrower, which is equivalent to spreading in k±. 



To summarize, we have proposed an efficient numerical method to solve Eq. 
via the Mellin transform. The method allows to study formal aspects of the equa- 
tions, e.g. the asymptotic forms of the evolution-generated UPDs at large 6, or at 
x — > and x — > 1. At large b the evolution-generated 6-dependent UPDs exhibit 
power-law fall-off, with the magnitude of the exponents growing with the probing 
scale [2]. The fall-off is steeper for the gluons than for the quarks. At x — > we have 
found generalizations of the DLL A behavior. We have also shown that for large b 
the solution for the valence UPD of the pion grows linearly with x for not too large 
x, and the slope decreases with b as a power law. At x — > 1 the evolution-generated 
6-dependent UPDs approach the integrated distributions as (1 — x) 2 . We find the 
spreading of the k± distributions with the probing scale Q, with the effect strongest 
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Figure 4. The rms transverse momenta of UPDs of the pion for x = 0.01, 0.1, and 0.5, plotted as 
functions of the renormalization scale Q 2 . Solid lines - gluons, dashed lines - non-singlet quarks, 
dotted lines - singlet quarks. 



for gluons and increasing with decreasing x. We have also shown that the widths 
(&5_)^ vo1 in all channels i increase at large Q 2 as Q 2 a(Q 2 ). 
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